Slow polaritons with orbital angular momentum in atomic gases 



J. Ruseckas, A. Mekys and G. Juzeliunas 
Institute of Theoretical Physics and Astronomy, Vilnius University, 
A. Gostauto 12, LT-01108 Vilnius, Lithuania 
(Dated: January 12, 2013) 

Polariton formalism is applied for studying the propagation of a probe field of light in a cloud 
of cold atoms influenced by two control laser beams of larger intensity. The laser beams couple 
resonantly three hyperfine atomic ground states to a common excited state thus forming a tripod 
configuration of the atomic energy levels involved. The first control beam can have an optical vortex 
with the intensity of the beam going to zero at the vortex core. The second control beam without 
a vortex ensures the loseless (adiabatic) propagation of the probe beam at a vortex core of the first 
control laser. We investigate the storage of the probe pulse into atomic coherences by switching 
off the control beams, as well as its subsequent retrieval by switching the control beams on. The 
optical vortex is transferred from the control to the probe fields during the storage or retrieval of 
the probe field. We analyze conditions for the vortex to be transferred efficiently to the regenerated 
probe beam and discuss possibilities of experimental implementation of the proposed scheme using 
atoms like rubidium or sodium. 
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FIG. 1. Probe and control laser beams acting on atoms characterized by two hyperfine ground states |1) and |2), as well as an 
excited state |0) to form a three-level scheme of the A type. Atoms are initially in the ground state |1) . Stimulated exchange 
of photons between the probe {£) and contol (f2c2) laser fields creates a superposition of the hyperfine atomic ground states |1) 
and 1 2) making the medium transparent for the resonant probe pulse. 
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FIG. 2. (a) Tripod scheme of the atom-light coupling involving a probe beams {£) and two control beams {Q,c2 and Qcs)- 
The three beams induce transitions between the atomic excited state |0) and three ground states |2) and |3). (a), (b) A 
control beam with the Rabi frequency f2c2 can have an optical vortex. Application of an additional control laser beam without 
the OAM {Q,c3 7^ 0) makes it possible to avoid losses in the propagation of the probe beam at the vortex core where the 
amplitude Qc2 goes to zero, (c) A possible experimental realization of the tripod setup for atoms like Sodium '7] or Rubidium 
[§] containing the hyperfine ground states with F = 1 and F = 2. 



I. INTRODUCTION 



During the last several years there has been a great deal of interest in slow [l|-|5|, stored |6l-ll8l] and stationary 
[l9l - l25j light. Light can be slowed down by seven orders of magnitude to velocities of several of tens of meters per 
second [l| due to the electromagnetically induced transparency (EIT) [26l - l30| . The EIT makes a resonant and opaque 
medium transparent for a probe beam by applying a control laser beam of larger intensity. The probe beam couples 
resonantly the ground and excited atomic states, whereas the control beam couples the same excited state to another 
unpopulated atomic ground state. This makes a A configuration of the atomic energy levels involved, as depicted 
in Fig. [T] The optical transitions induced by both laser beams interfere destructively preventing population of the 
excited atomic state. As a result, a weak pulse of probe light travels slowly and with little losses in a resonant medium 
due to the application of the control laser beam. 

The EIT was shown not only to slow down dramatically light pulses [I'Js'] , but also to store them 0, H, [l^ [TB - fTsj 
in atomic gases. The storage and release of a probe pulse has been accomplished 0, H, [13, ITsl - flsl by switching off 
and on the control laser 6] . The coherent control of the propagation of quantum light pulses can lead to a number 
of applications, such as generation of non-classical states in atomic ensembles and reversible quantum memories for 
slow light 6, 9, 10, 13, 29-33]. On the other hand, propagation of slow light through moving media 34-41] may be 
used for the light memories and rotational sensing devices. 

The orbital angular momentum (OAM) 'Jt', 'SsJ provides a new element to the slow light giving additional possi- 
bilities in manipulation of the optical information during the storage and retrieval of the slow light. The previous 
studies have concentrated on situations where the probe beam contains an OAM [4l| - |45| . In the present paper we 
consider another scenario in which it is a control laser beam which can carry an optical vortex. The intensity of such a 
control beam goes to zero at the vortex core leading to the absorption losses of the probe beam in this area. To avoid 
the losses we suggest to use an additional control laser without an optical vortex, so that the total intensity of the 
control lasers is non-zero at the vortex core of the first control laser. The probe and both control laser fields induce 
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transitions between the atomic energy levels in a tripod configuration of the light-atom coupling [49l - l60l | as depicted 
in Fig. ^!^. We show that the regenerated slow light can acquire the 0AM if one of the control beams contains it. 
We explore conditions for the optical vortex of the control beam to be transferred efficiently to the regenerated probe 
beam. 

The tripod scheme can be realised for atoms like Sodium ^] or Rubidium Q containing two hyperfine ground levels 
with F = 1 and = 2, as depicted in Fig. These atoms have been employed in the original experiments on the 
storage of slow light based on a simpler A setup 0, S|- the present situation |1) and |3) correspond to the magnetic 
sublevels (with mp = 1 and mp = —1) of the F = 1 hyperfine ground level, whereas the state |3) represents the 
hyperfine ground state with F = 2 and mp = 1- The probe beam is to be polarized, whereas both control beams 
are to be a~ polarized to make a tripod setup. Such a scheme can be produced by adding an extra circularly polarized 
laser beam fij-s as compared to the experiment by Liu et al 0] on the light storage in the sodium gases using the A 
scheme. Thus it is feasible to implement the suggested experiment on the transfer of optical vortex from the control 
to the probe fields using the tripod setup. 

The storage and retrieval of slow light is analyzed using the polariton formalism. The starting point is a set of the 
atomic equations together with the equation for the probe field. Subsequently we obtain two coupled equations for 
dark-state polaritons representing the slow light in the atomic medium. We provide conditions when the polaritons 
are decoupled. An advantage of polariton formalism is a simplicity of the relationship between the polariton field and 
the regenerated electric field, a feature which is missing in the direct analysis of the probe beam propagation fo^. 
Furthermore the equation for the polariton has a usual form of matter wave equation which describes the atomic 
evolution when the control fields are off. 



II. INITIAL EQUATIONS 

We will deal with an ensemble of atoms characterized by three hyperfine ground states |2), and |3), as well as 
an electronic excited state |0). The atomic internal and center of mass dynamics is represented by a four component 
field ^(r). Its components 'i>i{r,t), ^2(1", i), ^3(1", i), and 5'o(r,t) describe the atomic center of mass motion in the 
corresponding internal states |2), |3) and |0). In the semiclassical (mean field) approach, \['j(r,t) defines the 
probability amplitude to find an atom positioned at r in the j-th internal state, with j = 0,1,2,3. In the fully 
quantum approach ^!j(r,t) is the corresponding field operator. 

Three beams of light act on the atoms in a tripod configuration of the atom- light coupling [49l - l60j . Two strong 
classical control lasers induce transitions |2) — >■ |0) and |3) |0), whereas a weaker probe field drives a transition 
|1) — |0), as shown in Fig. [5] The former control lasers are characterised by the Rabi frequencies ilc2 and fics to be 
treated as incident variables. The latter probe beam is a dynamical quantity described by the electric field strength 

E(r,t)=e./^f(r,i)e~'"*+H.c., (1) 

where w = cfc is the central frequency of the probe photons, k = zfc is the wave vector, and e_Lz is the unit polarization 
vector. The probe field can be treated either as a classical variable or as a quantum operator. The dimensions of the 
electric field amplitude £ are chosen such that its squared modulus represents a number density of probe photons. 

The probe field is considered to be quasi-monochromatic, and its amplitude £ = £{r,t) changes little over the 
optical cycle. The slowly (in time) varying amplitude of the probe field obeys the following equation: 



where the parameter g = uj jle^h characterizes the strength of coupling of the probe field with the atoms, [l being 
the dipole moment of the atomic transition |1) — ?> |0). The quantities on the r.h.s. of Eq. ^ $0 and $* represent 
slowly (in time) varying atomic fields. The asterisk in <i>^ refers either to the complex conjugation of a classical field 
or the Hermitian conjugation of a quantum field. The slowly varying atomic fields $j (j — 1,2,3,4) are related to 
the original ones as: $1 = *ie*"i*, $2 = ^2e'("i+"""=2)*, $3 = ^ = 5'oe'('^i+'^)*, where fiwi is the 

energy of the atomic ground state 1, whereas Wc2 and Wcs are the frequencies of the control fields. 
The atomic equations of motion read 

A'$i = Fi(r)$i-%r$o, (3) 
k^Q = h{ujni ~ i7)$o + Vo(r)$o - ^^^2^2 - ?if^c3*3 - fig£^i , (4) 

A:$2 = /lW21«'2 + V'2(r)$2 - ftf^*2*0, (5) 
k^3 = huJ3i^3 + V3 (r)$3 - ftf^:3$0 , (6) 
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with 

where CJ21 = W2 ~ "^i + "^02 ^ a-nd W31 = W3 — + ui^'i — ^ are the frequencies of the electronic detuning from the 
two-photon resonances, wqi = ujo — uji ~ \s the frequency of the electronic detuning from the one-photon resonance, 
and 7 is the decay rate of the excited electronic level. Note that the inclusion of the non-zero decay rates should 
be generally accompanied by introduction of the noise operator in the equations of motion .281 . Yet in the present 
situation one can disregards the latter noise: we are working in the linear regime with respect to the probe field, 
so the population of the excited state is small enough. Here also m is the atomic mass and V^ (r) is the trapping 
potential for an atom in the internal state 7 (7 = 1,2,3,0). The terms containing atomic mass m are important for 
the description of the light-dragging effects [3J-[4l| . 

In Eqs. (121)~® the coupling of atoms with the probe and control fields has been written using the rotating wave 
approximation. Therefore, the last term in Eq. ([3]) has a negative frequency part of the probe field {£*), whereas 
the last term in Eq. (j4]) has a positive frequency part {£). Similarly Eq. (j4]) contains Rabi frequencies f2c2 and fics, 
whereas Eqs. ^ and ^ contain their complex conjugated counterparts 17*2 and $1*3. 

The equation of motion ([3]) for $1 does not explicitly accommodate collisions between the ground-state atoms. If 
the atoms in the internal ground-state 1 form a Bose- Einstein condensate (BEC), the coUisional effects can be included 
replacing Vi(r) by Vi(r) -I- guj^ip in Eq. ([3]), where 511 — 'inh? an / m and an is the scattering length between the 
condensate atoms in the internal state 1. This yields a mean- field equation for the condensate wave function $1. 

Initially the atoms populate the ground level 1. We are interested in the linear regime where the modulus of Rabi 
frequency of the probe field J7p = g£ is much smaller than the total Rabi frequency of the control beams fie given by 
Eq. © below. Consequently one can neglect the last term in Eq. ^ that causes depletion of the ground level 1. This 
provides a closed equation for the ground state dynamics: K^i = Vi{y)^i. If the atoms in the internal ground-state 
1 form a BEC, its wave function $1 = -^/n exp(iS'i) represents an incident variable determining the atomic density n 
and the condensate phase Si. The latter phase will not play an important role in our subsequent analysis, since we 
are not interested in the influence of the condensate dynamics on the propagation of slow light. The phase will be 
taken to be zero [Si — 0) when dealing with the storage and retrieval of slow light in the Section ITVl 

III. DARK- AND BRIGHT-STATE POLARITONS 

When the probe photons enter the atomic media, they are converted into composite quasiparticles of the radiation 
and atomic excitations known as polaritons. Let us first introduce the bright-state polariton $3: 



where 



is the total Rabi frequency. 



are dimensionless parameters, and 



*i3 =Cc(ec2$2+Cc3$3)+Clf, (8) 

= ^/Pc2? + Pe^\^' (9) 

U=^c2l^c. U = ^czl^c, (10) 

Ci=5'i>i/2, Cc = r!,/S (11) 



S = ^^11 + g^n . (12) 



The polariton $5 represents a specific superposition of the atomic and the probe fields featured in the equation 
of motion (U) for atoms in the excited electronic state. The latter equation (U) can be rewritten in terms of the 
bright-state polariton 

i^<I>o = h(ujm - i7)$o + yo(r)-I>o - K~<^b ■ (13) 

In this way the bright-state polariton is responsible for the light-induced atomic transitions to the excited state. 

The two dark-state polaritons are defined as superpositions of the atomic coherences and the probe photons or- 
thogonal to the bright-state polariton ^bi- 

<^Dl=Cc£-CliU^2+U^3), (14) 
$732 = C3*2 - Cc2^3 ■ (15) 
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It is to be noted that only the first dark-state polariton of the tripod scheme contains the electric probe field 
component and thus has a non-zero radiative group velocity. The incoming light is converted exclusively into this 
polariton when it enters the medium. The second dark-state polariton $£)2 does not have any contribution by the 
probe photons and is thus characterised by a zero radiative group velocity. It corresponds to the dark state of the A 
system consisting of the levels 2, 3, and 0. The combination ^£2^*2 + 'Cc3$3 featured in Eqs. (|8]) and (fT4)) . represents 
the bright state of such a A system. In this way, only the first polariton experiences the radiative motion, the second 
one being trapped in the atomic medium. 

The "bare" atomic and probe fields can be cast in terms of the dark and bright polaritons of the tripod system as: 



^"2 = C2(Cc$B - Ci^-ci) +?c34'r>2 



(16) 
(17) 
(18) 



To obtain the equation for the dark-state polaritons one needs to take the time derivative of Eqs. (ITi|) - ((T5|) and make 
use of the equations of motion ([2]), ([S])-® and (IT3|) . 

Suppose the control and probe beams are tuned close to the two-photon resonance. Application of such beams 
cause EIT in which the transitions |1) — |0), |2) — |0), and |3) — |0) interfere destructively preventing population 
of the excited state 0. As a result, the atom-light system is driven to the dark states, and the bright state polariton 
$B (featured in the equation of motion ([T^ for the excited state atoms) is weakly populated: $b ~ 0. Neglecting 
the contribution due to the bright-state polariton $b (adiabatic approximation), one obtains the equations for the 
dark-state polaritons and ^d2- Introducing a column $ = {^01,^02)^ , it is convenient to represent these 
equations in a matrix form: 
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l/mm 
1/m 



ihJ -W + U 



(19) 



where the 2x2 matrices J and U are defined in the Appendix!^ The former J represents a complex vector potential, 
U being a complex scalar potential. Even though the potentials are complex, the equation of motion ()19[) is Hermitian 
and thus it preserves the norm of the column $. Here also 



,2 



mDi= I ^ICcr + -|Ci|' 



(20) 



is the effective mass of the first dark-state polariton. The mass to_di exhibits position- and time-dependence through 
its dependence on the Rabi frequencies of the control fields and also on the atomic density. The second polaritons 
does not have a radiative component, so its effective mass coincides with the atomic mass m in Eq. (jl9p . 
The effective mass of the first polariton can be represented as: 



rriDi 



1 



TOrad m fl'^ + g'^n 



(21) 



where 



TOrad = =m (22) 

CWrad Wrad 



and 



+ g'^n 



^^ad = ^2 , (23) 



are, respectively, the radiative "mass " and the radiative group velocity of the first polariton, Vi-cc — huj/ nic being the 
atomic recoil velocity. In the slow light regime where Vtj. ^ g^n, the latter v^ad ~ (?Vij./g^n is much smaller than 
the vacuum speed of light: Wrad ^ c. The radiative velocity Wrad can be of the order of lOm/s for the slow light in 
atomic gases [1|. This greatly exceeds the typical velocities associated with the centre of mass motion of cold atoms. 
For instance, the atomic recoil velocity is typically of the order of 1 cm/s. Thus the second term can be neglected in 
Eq. ([HI), giving mm « mrad- 
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A. Co-propagating probe and control beams 

Suppose that the control beams propagate along z axis with krji ~ fcc3 = ^c- 

17e2 = ^^'^ , r!c3 - Kse'"'' ■ (24) 

For paraxial control beams the amplitudes f2^2 f-i^d ^c3 depend weakly on the propagation direction z. It is convenient 
to represent the dark-state polaritons as: 

$Di(r,t) = $'^i(r,i)e^'=% (25) 
$D2(r,t) = $'^2(r,i)e-^'=\ (26) 

where the amplitudes $^j^(r,t) and ^'^jI'^jO depend slowly on the propagation direction z in the paraxial case. 
Introducing a column <&' = {^'dii ^'02)^ ' ^'i- CS]), provides the following equation for the slowly varying amplitudes: 



ih 



d ^( vg, \ d 



y ( '"^1 ? ) V2 + thJ' • V + [/' 



at \ J dz 

where the 2x2 matrices J' and U' are presented in the Appendix [BJ Here 

h 



$' , (27) 



Wgl =«rad + -(fc-A:c)|Cl| (28) 

m 

is the group velocity of the first dark-state polariton. It comprises the radiative group velocity and the velocity of the 
two photon recoil. The latter term can be neglected giving Vgi ~ Wrad- 

B. Decoupled dark-state polaritons 

Let us analyse the terms which couple both dark polaritons in the equation of motion ()27p . The term with time 
derivatives in the non-diagonal elements of the matrix U' is proportional to 

If both control pulses depend on time in the same manner, i.e. r2c2 = fi^) ^^'^ "c3 = ^^c3 fi^)y above term 
is zero. Thus the coupling between the two dark-state polaritons can be avoided by switching both control pulses off 
and on in the same way, so that both of them exhibit the same temporal behaviour. 

Let us next estimate non-diagonal terms which contain the spatial derivatives of the control pulses in the equation 
of motion p7p and hence couple both dark polaritons. Such non-diagonal matrix elements are of the order of the 
atomic recoil energy fierce = ^^fc^/(2m) and thus can be neglected if the characteristic interaction time between the 
two dark-state polaritons Tpuisc = V^gi is small compared with the reciprocal recoil frequency: Wrcc'T'puiso ^ 1, where I 
is the length of the probe pulse in the medium. The latter condition can be easily fuUfillcd for typical slow light pulses 
whose durations are of the order of a microsecond p| and thus are much smaller than the reciprocal recoil frequencies. 
Consequently the polaritons $di and ^d2 are decoupled and equations for them can be solved separately. 

We are interested in the equation for the first dark polariton. Such a polariton contains the radiative contribution 
and thus describes propagation of the probe pulse of light in the medium. Neglecting the coupling with the second 
polariton, Eq. (1271) yields a closed equation for the paraxial propagation of the first polariton along the z direction: 

with Vgi ~ Wrad- Duc to the finite lifetime of the excited atomic state 7~^, the first polariton will experience radiative 
losses which are not included in the propagation equations (|27p and ((30|) . Let us now estimate the losses. The 
polariton lifetime is determined by the rate of the excited state decay and the total Rabi frequency of the control 
lasers ilc [§]: Tpoi = 7~^(Slc/Aa;)^, where Auj is a detuning from the two-photon resonance. One of the reasons for the 
appearance of the two-photon detunning is the finite duration of the probe pulse, Auj — Tp^^j^^,. To avoid the losses, a 
time the polariton tranverses the sample should be smaller than the polariton lifetime: i/wrad <SC ''"poi, with L being 
the length of the atomic cloud. This means the total Rabi frequency flc should be large enough, 

i<Wrad7"'f^c^pulsc- (31) 
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Note also that in the slow light regime, the probe radiation makes a tiny contribution to the polariton which is com- 
posed predominantly of the atomic excitations (atomic coherences). In fact, the velocity ratio Vrad/c ^ 1 represents 
a fraction of the radiation component in the polariton 0, . Thus Eq. (|30|) effectively describes propagation of the 
atomic coherences along the z axis at the velocity Wiad ^ c appearing due to the small radiative component. 



IV. STORAGE AND RELEASE THE SLOW LIGHT: GENERAL 



A. Storage of slow light 

Let us first consider the storage of the slow light. The probe beam enters the atomic medium at z — zq. The 
medium is illuminated by two control beams characterized by Rabi frequencies Q^^2 ^^'^ ^c3^ i where the index (s) 
refers to the storing stage of light. Initially the Rabi frequencies of the control beams (and hence the group velocity 
Vgi = Vg^i) are time- independent. Neglecting the diffraction effects, one can thus write 

£(^)(t,z)=£(^)(r,zo), T = t- r (l/vl'Adz'. (32) 

J Zo 

At the boundary the probe beam is converted into a dark-state polariton ^^^\{t) propagating at the group velocity 

(s) 

Vgi ^ c in the medium. Since the atomic population is created exclusively by the incident probe light, only the first 
dark-state polariton is populated, giving 

$(^) ^^(.)/^(.)^ $(^)^0, (33) 

where the temporal and spatial dependence of the first polariton are kept implicit. In writing the last relationship 
the use has been made of Eq. (|18p relating £ to and $5, together with the adiabatic approximation implying 

that ~ 0. For slow light the parameter Cc*^' « ^ v^]^ jc featured in Eq. ((33)) is much smaller than the unity. That's 

(s) . 

why the dark-state polariton contains only a tiny contribution by the electric field. 

The equations (fT6 |) ~([T7 |) together with the condition $^''(i) = 'I'^2(^) = ^ provide the atomic fields (atomic 
coherences) associated with the first polariton: 

- -df cf , <i>i^^ = -c^f d^^*^! . (34) 

At a certain time t = t^^^ the whole probe pulse enters the atomic medium and is contained in it. To store the slow 
polariton, both control fields are switched off at t = i^^^ in such a way that the Rabi frequency ratio r2|.2 /^ia = Cc2VCc3 
remains constant, whereas cj^"* 1- This gives the following atomic fields (atomic coherences) at the storing time: 

^ -e(r$(^) (t(^)) , $(^)(i(^)) ^ -^(f $(^l(i(^)) . (35) 

The stored atomic coherences no longer have the radiative group velocity and thus are trapped in the medium. The 
retrieval of these coherences is accomplished at a later time t = t^^^ . 



B. Regeneration of slow light 

To restore the polariton propagation, both control fields are switched on again at t = i*^''-' in such a way that their 
the ratio ^^2/^03 — ^c2 /^c3 constant. The difference between the storage and the retrieval times should not 
be too large, so that the atomic coherences given by Eq. ([55]) are preserved up to the retrieval time. In the initial 
experiment Q the light was stored up to a millisecond, yet the storage duration was increased up to a second recently 

(r) (r) (s) (s) 

If the relative Rabi frequencies ^^'^ C3 differ from the original ones ^^2 ^^'^ CcS > both dark-state polaritons 
are regenerated. Using Eqs. (fHl) and (IT5t . the dark state polaritons regenerated from the atomic coherences (I55t read 

(r) 

at the beginning of the release of light where Q sa 1: 

= iiilkit + ^^S^SKit^'^) , (36) 
= - ^iT^h^^olit^'^) ■ (37) 
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The electric probe field reappears due to the first dark-state polariton containing a non-zero electric field contribution: 

£:('-)(i) = CrH0<f^1W- (38) 

Substitution of Eq. ((36)) into Eq. (|38p and using Eq. (|33p , one can relate the regenerated electric field to the initial 
one as 



(r) 



(39) 



If both the storing and the retrieval takes place in the slow light regime, flc <C gy/n and flc <C g^/n , the above 
equation simplifies to 



(-,(r)^(s)* 
"c2 "c2 



^(r)f^(s)* 



Propagation of the regenerated polariton is governed by Eq. ([50)1 in the paraxial case. The polariton 

(r) 

propagates at the velocity Vgi = Vgi and might experience diffraction effects due to the second order transverse 
derivatives featured in Eq. ((30)) . On the other hand, the second polariton <f>£)2 is not coupled to the light fields and 
hence remains trapped (frozen) in the medium. 



(40) 



(r) 



V. STORAGE AND RETRIEVAL OF SLOW LIGHT: SPECIFIC SITUATIONS 

A. Restored control beams with the same spatial behaviour 

Let us first analyze the simplest situation where the Rabi frequencies of the restored control beams are proportional 
to the corresponding original ones with the same proportionality constant b: 

= , = bn^] (41) 

and hence ^^2'' — ^^'^ ^cs"* ~ ^(^3 ■ Under these conditions, Eqs. (|55)) and (|57)) together with (fTU)) provide the 
following amplitudes of the regenerated dark-state polaritons: 

$«(i«) = $(^l(i(^)), $W^o. (42) 

Thus the second polariton is not populated ($^2 ~ whereas the first regenerated dark-state polariton coincides 
with the original one. The corresponding regenerated electric field 

£;('■) = 6f (43) 

is proportional to the original one and thus does not acquire the phase singularity of the control beam r2c2 (if any) . In 
such a situation the vortex can not be transferred from the control to the regenerated probe beam. In the following 
Subsections we will analyse the vortex transfer from the control beam ilc2 to the regenerated probe beam in the case 
where the condition ()41|) no longer holds. Such a vortex transfer is accompanied with some population of the second 
polariton. 

It is noteworthy that the regenerated electric field given by Eq. (|15)) is increased (decreased) if the ratio of 

the total Rabi frequencies 6 = Jli'^-'/ric"'' is larger (smaller) than the unity. On the other hand the group velocity is 
increased for 6 > 1 and decreased for 6 < 1. This leads to the compression (for 6 < 1) or decompression for (6 > 1) of 
the regenerated probe pulse as compared to the stored one, a feature known from the light storage and retrieval in 
the A system Q. Note also that the total number of the regenerated photons is the same as that in the input beam. 
This is because the second polariton is not populated <I>£)2 = 0, so no atomic coherence remains frozen in the medium. 



B. Transfer of optical vortex at the retrieval of the probe beam 

(s) 

Suppose that only one control field is used during the storage phase of the probe light, i.e. fl).^ = and hence 
\^c^\ — 1- This means the storage stage involves a A scheme depicted in Fig. [T] In such a setup, the control beam 
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n^^2 can not carry an 0AM: Otherwise there would be non-adiabatic losses of the probe beam at the vortex core of 
the control beam. On the other hand, the retrieval of the probe beam is accomplished using a tripod system in which 
generally both fl^^^ ^^'^ ^^c3 ^^e non-zero. Under these conditions, Eqs. and ([57)) provide the following results 
for the regenerated polaritons: 

$S(ir^)=e^2^e^f'i>^i(i?), (44) 

^^DM^) = -iiTdt^^D\it^;^)- (45) 

The second polariton given by Eq. (j45p does not have a radiative component and is trapped in the medium. The 
electric field of the probe beam is regenerated exclusively due to the first polariton and reads using Eq. (j40|) 

(r) 

= (46) 

The equations (H^ or (|46p represent the initial condition for the subsequent propagation of the polariton in the 
medium. Such a polariton will propagate along the z axis with the group velocity Wrad- Its transverse profile will 
change due to the diffraction effects represented by the second order spatial derivatives in Eq. ([50)) . 

If the control beam ^lc2 carries an optical vortex at the retrieval stage, $7^2'* ~ e*^"^, the regenerated electric field 
^ gt^v acquires the same phase as one can see from Eq. (|46| . This means that the restored control beam transfers 
its optical vortex to the regenerated electric field £^'^\ If the initial control and probe fields have the same transverse 
dependence, the transverse profile of the regenerated electric field will mimic that of the control field f '^''^ ~ ri^j"* ~ c*^*^. 

As an illustration, let us take the restoring control laser f2^2^ to be the first order Laguerre-Gaussian (LG) beam: 

r2^2^ = ^pe*'^exp(-pVf^r)> where p = p/X is a dimensionless cylindrical radius, A — 2TT/k being the optical wave- 
length. On the other hand, the control beam is assumed to be the zero-order LG beam during the storage stage 
involving a A system: D,^^ = a~^Aexp{—p^/a'^) , where a determines a relative amplitude of the control fields fi^j' 

(s) 

and n),2 , Or and CTs being their dimensionless widths. This provides the following regenerated probe field 

= ape'^exp {-(? {a^^ - cj-^)] £^''>{&'>) . (47) 

It is noteworthy that the Rabi frequency of the additional laser fi^g^ does not enter the above equations (| ^ - (|T7)) for 
the regenerated probe field. Yet the additional laser plays an important role to ensure the lossles propagation of the 
restored probe field in a vicinity of the vortex core where ~^ 0) ^ o^^^ can see from Eq. (|31|). 



C. Transfer of the optical vortex during the storage of slow light 



Consider next the opposite situation where both control fields are on during the storage phase, so the storage of 
the probe beam is carried out using a tripod scheme. On the other hand, a A scheme is employed during the retrieval 
of the probe beam where only one control field is on, i.e. rj^3^ = and hence |^^2^| =1. In that case Eqs. (1361) and 
([37| yield the following results for the regenerated polaritons: 



c2 S.c2 
';c2 Sc3 



(48) 
(49) 



Again the electric probe field is regenerated exclusively due to the first polariton and is given by using Eq. (PO)) 



sir) 



^(r)^(s)* 
"c2 "c2 



(50) 



The equations (|15|) or ([501) represent the initial conditions for the subsequent propagation of the regenerated polariton 
governed, in the paraxial case, by the equation of motion ()30|) . Such a polariton will propagate along the z axis with 
the group velocity Wrad; and its transverse profile will change due to the diffraction effects represented by the second 
order spatial derivatives in Eq. pop . On the other hand, the second polariton will be frozen in the medium 
(neglecting the atomic motion) and its spatial form is given by Eq. (j49p . 
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If the second control beam carries an optical vortex at the storing stage, fJ^j e^^'^ , the regenerated electric field 

£ ~ e"*^''' acquires an opposite vorticity, as one can see from Eqs. (|48p and (|50p . The additional control beam fi^.g' 
does not have a vortex and hence is non-zero at the center. This ensures the lossless (adiabatic) propagation of the 
probe beam during the storage phase. It is noteworthy that the transverse profile of the regenerated probe field differs 
now from that of the storing beam due to the denominator in Eq. (|50p . 

Suppose that the control lasers are the first and zero order LG beams at the storage stage: 

= A~pe^^ exp(-pV'^') : ^^^3^ = Mexp(-pVa2) , (51) 

where the parameter h determines the relative amplitude of the additional control laser. On the other hand, the control 
beam is assumed to be the zero-order LG beam at the retrieval stage involing the A scheme: — a^exp(— 
Thus one arrives at the following regenerated probe field containing the phase conjugated vortex 

^^^^ = -2^/5e-^^exp [-p2 - .7^)] E^^\t^^^) . (52) 

It is noteworthy that for the 5 < 1 the transverse profile of the regenerated beam can differ considerably from the the 
Laguerre-Gaussian shape. Decreasing 6 the transverse shape of the regenerated beam becomes narrower. This leads 
to a larger difraction in its subsequent propagation, as it will be explored in the following Subsection. 



D. Dynamics of the restored probe beams 

Let us suppose that the atomic cloud is small enough, so that the diffraction can be neglected during the propagation 
of the probe beam in the medium. Such a condition can be fulfilled readily for a typical cloud of cold atoms, the 
length of which normally does not exceed a third of the millimeter [J. After leaving the atomic cloud, the probe 
beam propagates in the free space according to Eq. ^ with 5 = 0. Since the probe field is quasimonochromatic, its 
amplitude f (r, t) changes little during an optical cycle. In the stationary case one arrives at the following propagation 
equation for the slowly varying amplitude of the probe field: 

In the previous two subsections we have considered two possible scenarios to regenerate the probe field. In the 
first (A-T) case the Lambda scheme is used for storing the probe field whereas the tripod setup is employed for the 
regeneration. In the second (T-A) case, the tripod scheme is used for storing the probe field whereas the Lambda 
setup is used for the regeneration. In what follows we shall explore the subsequent propagation of the probe beam. 
The regenerated fields given by Eqs. (|47l) and ([5^ represent the initial conditions for such a propagation. By taking 
the initial probe beam to be Gaussian E'^"^ — ^q*' exp(— p^/dp), the regenerated fields read for both cases 

= a£!i'p^-e-^l^ , = .^^e^^ p,-^. 1^^ , (54) 

where = (7^^ + cr^^ — crj^ determines the width of the regenerated probe field measured in optical wavelength 
A = 27r/fc. 

The equation ([53l) has been solved numerically for a = 10. Figure [3] shows the subsequent propagation of the 

(r) 

regenerated beam for the first case. Here the regenerated field £j^_rp represents the first order LG beam and is 

(r) (s) 

proportional to the relative intensity of the control field at the release and storage stages f2^2 / ^c2 • The subsequent 
propagation of the field qualitatively preserves the transverse profile and is accompanied by some the diffraction 
spreading. 

The next three figures (Fig. 0] a,b,c) illustrate the second case where the tripod scheme is used for the storing and 

(r) 

the A scheme is used for the retrieval of the probe beam. The transverse profile of the regenerated beam £rp_/^ is 

(r) 

determined by the relative intensity b of the additional laser beam fi^g . For 6 = 3 the profile is much narrower as 
compared to the first order LG beam with the same width cr, as one can see comparing Figs. [3] and |4^. Consequently 
the light beam spreads out much faster than that in the first case (see Fig. |3^). When b increases (6 = 10 in Fig. |3Jd 
and 6 = 30 in Fig. [4j;), the profile of the probe beam approaches the shape featured in the first case. Note that the 
increase in the relative intensity b of the additional control laser is accompanied by the decrease in the intensity of the 
regenerated probe beam. Thus the improvement in the quality of the regenerated beams comes at a price of reducing 
its intensity. 
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FIG. 3. Propagation of the regenerated probe field in the free space for a = 10 and a = 1. The A scheme is used for storage 
and the tripod system for retrieval of the probe field. 



Polariton formalism has been applied for studying the propagation of a probe field of light in a cloud of cold atoms 
influenced by two additional control laser beams of larger intensity. The probe and control beams couple resonantly 
three hyperfine ground states to a common excited state in a tripod configuration of the atomic energy levels. The 
first control beam can have an optical vortex. Application of another control beam without a vortex ensures the 
loseless (adiabatic) propagation of the probe beam at the vortex core where the intensity of the first control beam 
goes to zero. The adiabatic propagation of the probe beam is obtained when the total intensity of the control lasers 
is sufficiently large at the vortex core. 

We have started with a set of atomic equations coupled with the equation for the probe beam, subsequently 
transforming them into two coupled equations for the dark-state polaritons. We have analysed conditions (related to 
the laser pulse durations and switching times) when the polaritons are decoupled and thus the problem reduces to 
a single equation for the polariton. An advantage of polaritonic analysis is a simplicity of the relationship between 
the polariton field and the regenerated electric field. Furthermore the equation for the polariton has a usual form of 
matter wave equation which describes also the atomic evolution when the control fields are off. 

The probe pulse is stored onto the atomic coherences and subsequently retrieved by switching off and on the control 
beams. As a result, the optical vortex can be transferred from the control to the probe fields during the storage or 
retrieval. Two scenarios have been analyzed in more details. The first case involves a A system for the storage and a 
tripod system for the retrieval. In such a situation the phase of vortex is transferred from the restoring control beam 
to the regenerated probe beam. In the second case the tripod system is used for the storage and the A system for the 
retrieval. The vortex phase is then transferred from the storing control beam to the regenerated probe beam in the 
phase conjugated form, so the probe beam acquires an opposite vorticity. The profile of the regenerated probe field is 
well preserved in the first case. On the other hand, in the second case the regenerated beam becomes narrower and 
thus experience larger difraction spreading. The width of the regenerated beam can be controlled by changing the 
intensity of the additional control beam during the storage phase. 

The tripod setup can be realized for atoms like Sodium Q or Rubidium containing two hyperfine ground levels 
with F = 1 and F = 2, as depicted in Fig. [2j;. The scheme can be produced by adding an extra circularly polarized 
laser beam f2c3 as compared to the experiment by Liu et al on the light storage in the gases using the A scheme. 
Thus it is feasible to implement the suggested experiment on the transfer of optical vortex from the control to the 
probe fields using the tripod setup. 



VI. CONCLUDING REMARKS 
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b) 










FIG. 4. Propagation of the regenerated electric field in the free space in the case of tripod system for storage and A system 
for retrieval by taking a — 10 and a — 1. The parameter b appearing in Eq. (I54|l is fe = 3, 10, 30 in the Figures a) b) and c), 
respectively. 
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Appendix A: Matrix elements in the equation for the dark polaritons 



The elements of the matrix J featured in Eq. (|T9|) . are 



Jii 



th ( — CcVCc 

ncu 



-crvci 



ICiPJb2: 



J22-*-(C3Vfc3+C2V^c2), 

m 

Ji2 = t-CUU^U - U^U) , 



(Al) 
(A2) 
(A3) 



where 



rB2 =*-(Cc2VC2+Cc3VC3) 



(A4) 
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The elements of the matrix U in Eq. p9)) . read 

/iV c2 ^ 1 



where 



Un = ( ^CcV^Cc + -CrV^Ci + ^ftC^VCl • Jb2 (A5) 

U22 - -|^(C3V'?c3 + C2^X2) (A6) 



U12 - -|^Cr(^c3V2Ce2 - U^'U) (A7) 
+ CrCc2Cc3(fi^32 + V^3(r)-1^2(r)) 

* f d d 

t^2i = -|^Ci (C2V'C3 - C3V'C2) + ^ft^vci • J2I (A8) 

+CiC2C3(^^32 + ^3(r)-V^2(r)) 

+*?^Cl fe3|C2-C2|C3 



t^B2 = -^i^c2^^C2 + ^czV'Cs) (A9) 

f d d 

+*MC2^ec2+C3^ec3 



Appendix B: Matrix elements in the paraxial equation for the dark polaritons 

Using Eq. ([M)) . the parameters ^c2 and ^c3 have the form ^c2 = ^'02^^^"^ ^ Cc3 = ■Ccs^*'^''^' where and slowly 
change with the distance z within the optical wavelength. The diagonal elements of the matrix J' entering Eq. <\27} . 
are given by 

J'n = fl^CcVCc + -CrVCi) + |CiPJb2 , (Bl) 



m 



where 



J22-*-(C3Vei3+C2Ve:2), (B2) 



JB2=*-(ec2VC2+e:3VC3)- (B3) 



The non-diagonal matrix elements of J' read 



J^^ = J2ie'('=+'=<=)^ = J'i*2 (B5) 



14 



The diagonal matrix elements of the complex scalar potential U' in Eq. (|27)) . are 

= f^-^-^'Cc + -CrV^i) + thQVCi ■ J'b2 (B6) 
2 \riu! m J 

C^22 = -|^(C3V^^^3 + C2^'C2) (B7) 



(fiL^21 + "t^2(r))|ee3r + (^^31 + V^irMc 

■ ( t d I* id 

( '?c3^Cc3 + ?c2^Cc2 ) > 



where 



C/b2 - -^{^2^'^ + e:3V^C3) (B8) 
+ (fiw21 + ^2(r))|^c2|' + (fiW31 + T/3(r))|Cc3|' 

■ ( I* d I I* 9 I 

I Cc2^Cc2 + Cc3^^c3 



Finally, the non-diagonal elements of the complex scalar potential U' are given by 

U[2 = -^Cie'^'^-'^'-iC3^X2 - C2^'C3) 

2m 



+Cre^^^=-'=)^ec2C3(fi^32 + ^3(r) - V2{V)) 



Kc2q^Sc3 Sc3 ^^^02 ; ' 

zm (,1 



+Cie»('=-'=^'^C2a(^32 + V3(r) - V2{V)) 
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